#######################
# World              ##
#######################
d:Ge/World/HLX   = 100 nm
d:Ge/World/HLY   = 100 nm
d:Ge/World/HLZ   = 100 nm
s:Ge/World/Color = "white"
s:Ge/World/Material = "G4_WATER"
b:Ge/World/Invisible = "True"

#######################
# NP
#######################
s:Ge/NP/Type               = "TsSphere"
s:Ge/NP/Parent             = "World"
s:Ge/NP/Material           = "G4_Au"
d:Ge/NP/TransX             = 0 nm
d:Ge/NP/TransY             = 0 nm
d:Ge/NP/TransZ             = 0 nm
d:Ge/NP/RotX               = 0 deg
d:Ge/NP/RotY               = 0 deg
d:Ge/NP/RotZ               = 0 deg
d:Ge/NP/RMin               = 0 nm
d:Ge/NP/RMax               = 5 nm
d:Ge/NP/SPhi               = 0 deg
d:Ge/NP/DPhi               = 360 deg
d:Ge/NP/STheta             = 0 deg
d:Ge/NP/DTheta             = 180 deg
s:Ge/NP/Color              = "yellow"

###############################################
#####          Physics                  #######
###############################################
s:Ph/ListName                                         = "Fisica"
sv:Ph/Fisica/Modules                                  = 1 "g4em-livermore"
s:Ph/Fisica/Type                                      = "Geant4_Modular"
d:Ph/Fisica/SetProductionCutLowerEdge = 100 eV
d:Ph/Fisica/CutForGamma           = 1 nm             # overrides CutForAllParticles for Gamma
d:Ph/Fisica/CutForElectron        = 1 nm             # overrides CutForAllParticles for Electron
b:Ph/Fisica/Fluorescence          = "True"           # Set to true to turn on Fluorescence
b:Ph/Fisica/Auger                 = "True"           # Set to true to turn on Auger
b:Ph/Fisica/AugerCascade          = "True"           # Set to true to turn on AugerCascade
b:Ph/Fisica/DeexcitationIgnoreCut = "True"           # Set to true to implement DeexcitationIgnoreCut
b:Ph/Fisica/PIXE                  = "True"           # Set to true to turn on PIXE

#######################
# Source
#######################
s:Ge/BeamPosition/Parent  = "World"
s:Ge/BeamPosition/Type    = "Group"
d:Ge/BeamPosition/TransX  = 0. nm
d:Ge/BeamPosition/TransY  = 0. nm
d:Ge/BeamPosition/TransZ  = -1 * Ge/NP/RMax nm 
d:Ge/BeamPosition/RotX    = 0 deg
d:Ge/BeamPosition/RotY    = 0. deg
d:Ge/BeamPosition/RotZ    = 180. deg

s:So/Phase_1/Type                            = "PhaseSpace"
s:So/Phase_1/PhaseSpaceFileName              = "PhaseSpace_1"
s:So/Phase_1/Component                       = "BeamPosition"
b:So/Phase_1/PhaseSpaceIncludeEmptyHistories = "True"

# Scaling factor = R_NP/R_Cylinder
u:So/Phase_1/PhaseSpaceScaleXPosBy =  5.26e-7
u:So/Phase_1/PhaseSpaceScaleYPosBy =  5.26e-7
u:So/Phase_1/PhaseSpaceScaleZPosBy =  5.26e-7

###############################################
#####          Outputs                  #######
###############################################

# Phase Space File
s:Sc/PhaseSpace/Quantity                    = "PhaseSpace"
b:Sc/PhaseSpace/OutputToConsole             = "True"
s:Sc/PhaseSpace/Surface                     = "NP/OuterCurvedSurface"
s:Sc/PhaseSpace/OutputType                  = "ASCII" # ASCII, Binary, Limited or ROOT
s:Sc/PhaseSpace/OutputFile                  = "Nanoparticle"
i:Sc/PhaseSpace/OutputBufferSize            = 1000
s:Sc/PhaseSpace/OnlyIncludeParticlesGoing  = "Out"
b:Sc/PhaseSpace/IncludeTOPASTime            = "False"
b:Sc/PhaseSpace/IncludeTimeOfFlight         = "False"
b:Sc/PhaseSpace/IncludeRunID                = "False"
b:Sc/PhaseSpace/IncludeEventID              = "False"
b:Sc/PhaseSpace/IncludeTrackID              = "False"
b:Sc/PhaseSpace/IncludeParentID             = "True"
b:Sc/PhaseSpace/IncludeCreatorProcess       = "True"
b:Sc/PhaseSpace/IncludeVertexInfo           = "False"
b:Sc/PhaseSpace/IncludeSeed                 = "False"
s:Sc/PhaseSpace/IfOutputFileAlreadyExists   = "Overwrite"
sv:Sc/PhaseSpace/OnlyIncludeParticlesNamed = 1 "e-"

# Deposited Energy inside the NP
s:Sc/EnergyDeposit/Quantity                         = "EnergyDeposit" # sum of energy deposits divided by mass
s:Sc/EnergyDeposit/Component                        = "NP"
b:Sc/EnergyDeposit/OutputToConsole                  = "False"
s:Sc/EnergyDeposit/IfOutputFileAlreadyExists        = "Overwrite"
sv:Sc/EnergyDeposit/Report         = 2 "Sum" "Standard_Deviation" 

# Deposited Energy inside the NP by secondary particles
s:Sc/EnergyDeposit_Sec/Quantity                         = "EnergyDeposit" # sum of energy deposits divided by mass
s:Sc/EnergyDeposit_Sec/Component                        = "NP"
b:Sc/EnergyDeposit_Sec/OutputToConsole                  = "False"
s:Sc/EnergyDeposit_Sec/IfOutputFileAlreadyExists        = "Overwrite"
sv:Sc/EnergyDeposit_Sec/Report         = 3 "Sum" "Standard_Deviation" "Count_In_Bin"
s:Sc/EnergyDeposit_Sec/OnlyIncludeParticlesOfGeneration = "Secondary"

#######################
# Graphics
#######################
b:Gr/ViewA/Active = "False"
s:Gr/ViewA/Type             = "OpenGL"
i:Gr/ViewA/WindowSizeX      = 1024
i:Gr/ViewA/WindowSizeY      = 768
b:Gr/ViewA/IncludeAxes      = "True"
d:Gr/ViewA/Theta            = 55 deg
d:Gr/ViewA/Phi              = 20 deg
s:Gr/ViewA/Projection       = "Perspective"
d:Gr/ViewA/PerspectiveAngle = 30 deg
u:Gr/ViewA/Zoom             = 1.
d:Gr/ViewA/AxesSize = 3 nm # size of axes
s:Gr/ViewA/BackgroundColor = "white" 
Ts/UseQt = "false"
b:Ts/PauseBeforeQuit = "False"
